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Abstract. A microscopic heterogeneous system under random influence 
is considered. The randomness enters the system at physical boundary of 
small scale obstacles as well as at the interior of the physical medium. This 
system is modeled by a stochastic partial differential equation defined on a 
domain perforated with small holes (obstacles or heterogeneities), together 
with random dynamical boundary conditions on the boundaries of these 
small holes. 

A homogenized macroscopic model for this microscopic heterogeneous 
stochastic system is derived. This homogenized effective model is a new 
stochastic partial differential equation defined on a unified domain with- 
out small holes, with static boundary condition only. In fact, the random 
dynamical boundary conditions are homogenized out, but the impact of 
random forces on the small holes' boundaries is quantified as an extra sto- 
chastic term in the homogenized stochastic partial differential equation. 
Moreover, the validity of the homogenized model is justified by showing 
that the solutions of the microscopic model converge to those of the effec- 
tive macroscopic model in probability distribution, as the size of small holes 
diminishes to zero. 
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1. Introduction 

Stochastic effects in the multiscale modeling of complex phenomena have 
drawn more and more attention in many areas such as material science |10j . 
climate dynamics [27], chemistry and biology [2TII51] . Stochastic partial differ- 
ential equations (SPDEs or stochastic PDEs) arise naturally as mathematical 
models for multiscale systems under random influences. The need to include 
stochastic effects in mathematical modeling of some realistic complex behaviors 
has become widely recognized in science and engineering. But implementing 
this approach poses some challenges both in mathematical theory and compu- 
tation [SI [TEl EH [27J EU H3] . The addition of stochastic terms to mathematical 
models has led to interesting new mathematical problems at the interface of 
dynamical systems, partial differential equations, scientific computing, and 
probability theory. 

Sometimes, noise affects a complex system not only inside the physical 
medium but also at the physical boundary. Such random boundary condi- 
tions arise in the modeling of, for example, the air-sea interactions on the 
ocean surface [32], heat transfer in a solid in contact with a fluid [31], chem- 
ical reactor theory [32], and colloid and interface chemistry [5B]. Random 
boundary conditions may be static or dynamical. The static boundary condi- 
tions, such as Dirichlet or Neumann boundary conditions, do not involve with 
time derivatives of the system state variables. On the contrary, the dynami- 
cal boundary conditions contain such time derivatives. Randomness in such 
boundary conditions are often due to various fluctuations. 

In this paper we consider a microscopic heterogeneous system, modeled by 
a SPDE with random dynamical boundary condition, in a medium which ex- 
hibits small-scale spatial heterogeneities or obstacles. One example of such 
microscopic systems of interest is composite materials containing microscopic 
holes (i.e., cavities), under the impact of random fluctuations in the domain 
and on the surface of the holes [281 ES] . A motivation for such a model is based 
on the consideration that the interaction between the atoms of the different 
compositions in a composite material causes the thermal noise when the scale 
of the heterogeneity scale is small. A similar consideration appears also in a 
microscopic stochastic lattice model [6] for a composite material. Here the 
microscopic structure is perturbed by random effect and the complicated in- 
teractions on the boundary of the holes is dynamically and randomly evolving. 
The heterogeneity scale is assumed to be much smaller than the macroscopic 
scale, i.e., we assume that the heterogeneities are evenly distributed. From a 
mathematical point of view, one can assume that microscopic heterogeneities 
(holes) are periodically placed in the media. This spatial periodicity with small 
period can be represented by a small positive parameter e (i.e., the period). In 
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fact we work on the spatial domain D € , obtained by removing S e , a collection 
of small holes of size e, periodically distributed in a fixed domain D. When 
taking e — > 0, the holes inside domain D are smaller and smaller and their 
numbers goes to oo. This signifies that the heterogeneities are finer and finer. 

In other words, we consider a spatially extended system with state variable 
u e , where stochastic effects are taken into account both in the model equation 
and in the boundary conditions, defined on a domain perforated with small 
scale holes. Specifically, we study a class of stochastic partial differential equa- 
tions driven by white noise on a perforated domain with random dynamical 
boundary conditions: 

du e (t,x) = Au e (t, x) + f(t, x, u t , Vw e ) dt + gi(t,x)dWi(t,x) 

in D e x (0,T), 
r diii (t x ) 

e 2 du t (t,x) = ^— ebu t (t,x) dt + eg2(t, x)dW2(t, x) 

on dS e x (0,T). 

This model will be described in more detail in the next section. 

The goal is to derive a homogenized effective equation, which is a new sto- 
chastic partial differential equation (see Theorems 15.11 16. 1\ 16.21 and 16.31) . for 
the above microscopic heterogenous system, by homogenization techniques in 
the sense of probability. Homogenization theory has been developed for deter- 
ministic systems, and compactness discussion for the solutions {u t } e in some 
function space is a key step in various homogenization approaches [12]. How- 
ever, due to the appearance of the stochastic terms in the above microscopic 
system considered in this paper, such compactness result does not hold for 
this stochastic system. Fortunately the compactness in the sense of proba- 
bility, that is, the tightness of the distributions for {u t }, still holds. So one 
appropriate way is to homogenize the stochastic system in the sense of prob- 
ability. It is shown that the solution u e of the microscopic or heterogeneous 
system converges to that of the macroscopic or homogenized system as e j in 
probability distribution. This means that the distribution of {u e } e weakly con- 
verges, in some appropriate space, to the distribution of a stochastic process 
which solves the macroscopic effective equation. 

It is interesting to note that, for the above system with random dynamical 
boundary conditions, the random force on the boundary of microscopic scale 
holes leads, in the homogenization limit, to a random force distributed all 
over the physical domain D, even when the model equation itself contains 
no stochastic influence in the domain; see Remark 15.21 in §5. We could also 
say that the impact of small scale random dynamical boundary conditions 
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is quantified or carried over to the homogenized model as an extra random 
forcing. Therefore, the homogenized effective model is a new stochastic partial 
differential equation, defined on a unified domain without holes. 

In the present paper, the two-scale convergence techniques are employed 
in our approach. Two-scale convergence method is an important method in 
homogenization theory which is a formal mathematic procedure for deriving 
macroscopic models from microscopic systems. Two-scale convergence method 
contains more information than the usual weak convergence method; see [2] or 
§U Moreover by use of the two-scale convergence, we do not need the extension 
operator as introduced in [T5] . 

Partial differential equations (PDEs) with dynamical boundary conditions 
have been studied recently in, for example, [U [201 1221 [23J, [25j H7] and reference 
therein. The parabolic SPDEs with noise in the static Neumann boundary 
conditions have also been considered in [T6J [T7J [36] . In [11] , the authors have 
studied well-posedness of the SPDEs with random dynamical boundary condi- 
tions. One of the present authors, with collaborators, has considered [T8J I5T] 
dynamical issues of SPDEs with random dynamical boundary conditions. 

The homogenization problem for the deterministic systems defined in per- 
forated domains or in other heterogeneous media has been investigated in, for 
example, [H |39l HOI HQ HH] for heat transfer in a composite material, [8], [13]. IT5] 
for the wave propagation in a composite material and [3^1 [38] for the fluid flow 
in a porous media. For a systematic introduction in homogenization in the 
deterministic context, see [12], [281 SSI In [S], the effective macroscopic 
dynamics of a deterministic partial differential equation with deterministic 
dynamical boundary condition on the microscopic heterogeneity boundary is 
studied. 

Recently there are also works on homogenization of partial differential equa- 
tions (PDEs) in the random context; see [23 EH EJ [28] for PDEs with 
random coefficients, and [3 EH [59j [28] for PDEs in randomly perforated 
domains. A basic assumption in these works is the ergodic hypotheses on 
the random coefficients, for the passing of the limit as e — > 0. Note that 
the microscopic models in these works are partial differential equations with 
random coefficients, so-called random partial differential equations (random 
PDEs)0 M, SD E2 EE E3], instead of stochastic PDEs — PDEs with noises 
- in the present paper; see also [52]. Another novelty in the present paper 
is that the microscopic system is under the influence of random dynamical 
boundary conditions. 

We first consider the linear system and then present results about nonlinear 
systems with special nonlinear terms. This paper is organized as follows. The 
problem formulation is stated in $2j Section [3] is devoted to basic properties of 



HOMOGENIZED DYNAMICS OF SPDES 



5 



the microscopic heterogeneous system, and some knowledge to be used in our 
approach is introduced in §U The homogenized effective macroscopic model 
for the linear system is derived in §0 In the last section, homogenized effective 
macroscopic models are obtained for three types of nonlinear systems. 



2. Problem formulation 

Let the physical medium D be an open bounded domain in R", n > 2, 
with smooth boundary dD, and let e > be a small parameter. Let Y = 
[0, h) x [0, 12) x • • • x [0, l n ) be a representative elementary cell in W 1 and S an 
open subset of Y with smooth boundary dS, such that S CY. The elementary 
cell Y and the small cavity or hole S inside it are used to model small scale 
obstacles or heterogeneities in a physical medium D. Write I = l 2 , • ■ ■ , l n ). 
Define eS = {ey : y G S}. Denote by S tt k the translated image of eS by kl, 
k G Z n , kl = (kih, k 2 h, • • • , k n l n ). And let S e be the set all the holes contained 
in D and D e = D\S e . Then D e is a periodically perforated domain with holes 
of the same size as period e. We remark that the holes are assumed to have no 
intersection with the boundary dD, which implies that 8D e = dD U dS t . See 
Fig. 1 for the case n — 2. This assumption is only needed to avoid technicalities 
and the results of our paper will remain valid without this assumption [3]. 




Fig. 1: Geometric setup in M. 2 

In the sequel we use the notations 

Y* = Y\S, *? = ^ 

with \Y\ and \Y*\ the Lebesgue measure of Y and Y* respectively. Denote by 
X the indicator function, which takes value 1 on Y* and value on Y \ Y*. In 
particular, let xa be the indicator function of A C M n . Also denote by v the 
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zero extension to the whole D for any function v defined on D e : 

~ _ j v on D t , 
V ~ \ on S e . 

Now for T > fixed final time, we consider the following Ito type nonau- 
tonomous stochastic partial differential equation defined on the perforated do- 
main D t in R n 



du t (t, x) 

e 2 du t (t, x) 

u e (t,x) 
u e (0,x) 



Au e (t, x) + f(t, x, u e , Vu e ) dt + gi(t, x)dWi(t, x) 

in D t x (0,T), 
du e (t, x) 



(2.1) 



on dS 6 x (0,T), 
on dD x (0,T), 
u (x) in D e , 



-ebu e (t,x) dt + eg 2 (t,x)dW 2 (t,x) (2.2) 



(2.3) 
(2.4) 

satisfies some prop- 



where b is a real constant, / : [0, T] x D x R x IR n — 
erty which will be described later and u e is the exterior unit normal vector 
on the boundary dS e , v G L 2 (dS e ) and u G L 2 (D). Moreover, Wi(t,x) and 
W2(t,x) are mutually independent L 2 (D) valued Wiener processes on a com- 
plete probability space (Q,^ 7 , P) with a canonical filtration (.Ft)t>o- Denote 
by Qi and Q2 the covariance operators of W\ and W 2 respectively. Here we 
assume that gi(t,x) G C(L 2 (D)), i = 1,2 and that there is a positive constant 
Ct independent of e such that 



\9i(t,-)\\: 



v — "\ — 



T- 



1,2, tG[0,T] 



(2.5) 



where {ej}'jL l are eigenvectors of operator —A on D with Dirichlet bound- 
ary condition and they form an orthonormal basis of L 2 (D). Here £(L 2 (D)) 
denotes the space of bounded linear operators on L 2 (D) and C 2 = £^ l (if) 
denotes the space of Hilbert-Schmidt operators related to the trace operator 
Qi [IE] . We also denote by E the expectation operator with respect to P. 

Let S be a Banach space and S' be the strong dual space of S. We recall the 
definitions and some properties of weak convergence and weak* convergence 



Definition 2.1. A sequence {s n } in S is said to converge weakly to s G S if 
Ws'eS', 

lim {s',s n ) S ',s = {s',s) S ',s 
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which is written as s n — ^ s weakly in S. Note that (s',s) denotes the value of 
the continuous linear functional s' at the point s. 

Lemma 2.2. (Eberlein-Shmulyan) Assume that S is reflexive and let {s n } 
be a bounded sequence in S. Then there exists a subsequence {s nk } and s G 
S such that s nk s weakly in S as k — > oo. If all the weak convergent 
subsequence of {s n } has the same limit s, then the whole sequence {s n } weakly 
converges to s. 

Definition 2.3. A sequence {s' n } in S' is said to converge weakly* to s' G S 1 
ifWseS, 

lim (s' n ,s) S ',s = (s',s) S ',s 

n— >oo 

which is written as s' n — s' weakly* in S' . 

Lemma 2.4. Assume that the dual space S' is reflexive and let {s' n } be a 
bounded sequence in S' . Then there exists a subsequence {s' nk } and s' G S' 
such that s' nk — ^ s' weakly* in S' as k — > oo. If all the weakly* conver- 
gent subsequence of {s' n } has the same limit s' , then the whole sequence {s ; n } 
waekly* converges to s' . 

In the following, for a fixed T > 0, we always denote by Ct a constant 
independent of e. And denote by D T the set [0,T] x D. 

3. Basic properties of the microscopic model 

In this section we will present some estimates for solutions of the microscopic 
model (12.1 1) , and then discuss the tightness of the distributions of the solution 
processes in some appropriate space. We focus our argument in the case of 
linear microscopic systems, where the term / is independent of u e and V« e and 
/(-, ■) G L 2 (0,T; L 2 (D)). Then we briefly extend this to the case of nonlinear 
microscopic systems with Lipschitz nonlinearities. 

Define by H^(D e ) the space of elements of H l (D e ) which vanish on dD. 
Denote by H~ l (D e ) the dual space of H}(D e ) with the usual norm and let 
7 e : if 1 (Z) e ) — > L 2 (dS e ) be the trace operator with respect to dS t which is 

continuous [32]. We also denote that H^(dS e ) = 7,(i/ 1 (A)) and let H~*(D,) 

i 

be the dual space of Hi (D t ). 

Introduce the following function spaces 

X] = {(u,v) G H](D e ) x Hl{dS e ) : v = e%u} 

and 

X° = {L 2 (A)xL^)} 
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with the usual product and norm. Define an operator B t on the space H*(D e ) 
as 

Oil 

B e u = — + ebu, u G H}(D e ). (3.1) 
on 

Now we define the operator A e on V(A e ) = {(u, v) G X] : (—Am, R t B e u) G 
A e }, where R e is the restriction to dS e , as 

A e z = (-Am, -R e B e u), z = (u,v) G V(A e ). (3.2) 
e 

Associated with the operator A e , we introduce the bilinear form on X\ 

a e (z,z) = / VuVudx + eb I 7 e (-u)7 e (w)<is (3.3) 

with z = («,m), z = (u,v) G X*. Notice that \le{u)\ 2 L 2i Qs ^ < CXSJIul 2 ^^, 
we see that there is M > 0, independent of e, such that 

a £ (z, z) < M\u\ H i {Dt) \u\ H i (Dt) 

and the following coercive property of a holds 

a e (z, z) > a\z\ 2 x i - f3\z e \ 2 x o, (3.4) 

for some constants a, (3 > which are also independent of e. Write the Co- 
semigroup generated by operator — A t as S e (t). 

Then the system (I2.1I) - (I2.4I) can be rewritten as the following abstract sto- 
chastic evolutionary equation 

dz £ (t, x) = [-A e z e (t, x) + F £ (t, x)]dt + G e (t, x)dW(t, x), z e (0) = z (3.5) 

where 

F e (t, x) = (f(t, x), 0)', G e (t, x)dW{t) = ( 9l (t, x)dW 1 (t, x),g 2 {t, x)dW 2 {t, x)f 
and z = (uo,vo). And the solution of (13.51) can be written in the mild sense 

Ze (t) = S e (t)z + [ S e (t - s)F e (s)ds + [ S e (t- s)G e (s)dW(s). (3.6) 
Jo Jo 

Moreover, the variational formulation is 

nii e (fdxdt + e 2 / / ii € (pdxdt + eb / u € (pdxdt = (3.7) 
_> E Jo JdS t Jo JdS t 

— j j VutVtpdxdt + fipdxdt + g\tpW\dxdt + e / / g2<pW2dxdt 

Jo Jd € Jo Jd c Jo Jd € Jo JdS e 

for (p(t,x) G C^{[0,T] x D e ). Here ' denotes |. 

For the well-posedness of system (13.51) we have the following result. 
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Theorem 3.1. (Global well-posedness of microscopic model) Assume 
that 1(2. 5\) holds for T > 0. Ifz = (uq, v ) is a (jF , B(X®)) -measurable random 
variable, then the system $3. 5\) has a unique mild solution z e E L 2 (Tl,C(0,T; 
X e °) fl L 2 (0, T; X*)) , which is also a weak solution in the following sense 

{ze{t)A)xl (3.8) 
= (zo,<P)x}+ / (-A e z e (s),<f)) x} ds+ / (F ej <f)) x ids+ / (G € dW,(j))x} 



for t E [0, T) and E X\. Moreover if zq is independent of W{t) with 
E|z |xo < °°j then 

E\z e {t)\ 2 x0 + E [ \z e {s)\ 2 x ids < (l + E\z \ x o)C T , forte[0,T] (3.9) 
Jo 



and 



E{ sup \z e (t)\ 2 x0 } < (l + E|z ||o+E / \z e (s)\ 2 xl ds)C T . (3.10) 
te[o,T] ' Jo 

Proof. By the assumption (12. 5p . we have 

\\G t (t,x)\\ 2 Q = \\gi(t,x)\\ 2 Ql + \\92{t,x)\\ 2 Qa < oo. 

i- 2 *-2 *-2 

Then the classical result P2] yields the local existence of z € . By applying the 
stochastic Fubini theorem |16| . it can be verified that the local mild solution 
is also a weak solution. 

Now we give the following a priori estimates which yields the existence of 
weak solution on [0, T] for any T > 0. 

Applying Ito formula to |z e |^. , we derive 

d\z e (t)\ 2 x o+2(A e z e ,z e ) x odt = 2(F e (t,x), z e ) x odt + 2(G £ (t,x)dW(t), z £ ) x o + 

\G e (t,x)\ 2 Q dt. (3.11) 

By the coercivity (13.41) of a e (-, •), integrating (13. lip with respect to t yields 

\z e (t)\ 2 x o + 2a / \z e (s)\ x ids 
Jo 

+ (2/3 + 1) / \z e (s)\ 2 x0 ds + 
Jo 

2 [ {G e {s)dW{s),z e (s)) x o+ [ \G e (s)\ 2 Q ds. 
Jo Jo 2 
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Taking expectation on both sides of the above inequality yields 

E\z e (t)\ 2 x0 + 2aE [ \z e (s)\ 2 x ids 
Jo 

< E\z \ 2 x0 + \F e \ 2 L2{0tT . x0) + (2/9 + 1) / E\z e (s)\ x0 ds + / \G e {s)\ 2 Q ds. 

Jo Jo 2 

Then the Gronwall lemma gives the estimate ( 13. 91) . Notice that, by Lemma 
7.2 in [TBI, 



E sup 

te[o,T] 



S^t — s)G e (s, x)ds 



T 

<CV / |G f fs)| 2 , ds. 



x - -'0 



Therefore by the assumption on / and (13. 6p we have the estimate (I3.10p . The 



proof is hence complete. □ 
By the above result and the definition of z e we have the following corollary. 

Corollary 3.2. Assume the conditions in Theorem \3.1i Then for t e [0, T], 

we have 

HK(t)\ 2 L2{Df) + e 2 \ 1 Mt)\l Hd s e )) + 

+ e 2 \^ e u e (t)\ 2 Hl/2{ds ^ds < (1 + E|z q |^o)Ct (3.12) 

J 

and 

E{ sup \u e (t)\h {Dc) + t 2 \lMt)\li {dSt )} < (! + Ekolxo)CT- (3.13) 

te[o,T] 



We recall a probability concept. Let z be a random variable taking values 
in a Banach space S, namely, z : — >• 5. Denote by C(z) the distribution (or 
law) of z. In fact, C(z) is a Borel probability measure on S defined as [16] 



C(z)(A) = F{u : z(u) e A}, 

for every event (i.e., a Borel set) A in the Borel cr— algebra B(S), which is the 
smallest a— algebra containing all open balls in S. 

As stated in £JU for the SPDE ( 12.1 ft we aim at deriving an effective equation 
in the sense of probability. A solution u 6 may be regarded as a random variable 
taking values in L 2 (0, T; L 2 (D € )). So for a solution w e of fl2.1D - (j2.4j) defined on 
[0,T], we focus on the behavior of distribution of u € in L 2 (0, T; L 2 (D e )) as 
e — > 0. For this purpose, the tightness [19J of distributions is necessary. Note 
that the function space changes with e, which is a difficulty for obtaining the 
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tightness of distributions. Thus we will treat {C(u e )} family of dis- 

tributions on L 2 (0,T; L 2 (D)) by extending u e to the whole domain D. Recall 
that the distribution for law ) of u f is defined as: 



C{u e ){A) = \ 
for Borel set A in L 2 (0, T; L 2 (D)). 



,to) E A} 



First we define the following spaces which will be used in our approach. 
For Banach space U and p > 1, define W 1,p (0,T; U) as the space of functions 
h G LP(0,T;U) such that 

dh p 



\h\ 



W 1 >p(0,T;£T) 



(It 



LP(0,T;U) 



< oo. 



And for any a G (0,1), define W a,p (0, T; U) as the space of function h G 
L p (0,T;U) such that 



\hf 



W a 'P(0,T;U) 



\LP(P,T;U) 



+ 



JO 



W)-h(s)\\ 
\f. - m\ 1 +°p 



-dsdt < oo. 



For p G (0, 1), we denote by C p (0, T; U) the space of functions h : [0, T] — > X 
that are Holder continuous with exponent p. 

Theorem 3.3. (Tightness of distributions) Assume that z = (uq,Vq) be 
a LFo, B(X®)) -measurable random variable which is independent ofW(t) with 
E\z \ 2 x0 < oo. Then for any T > 0, (C(u e )) e! the distributions of {u e ) e , is tight 
in L 2 (6, T; L 2 (D)) n C(0, T; H~\D)). 

Proof. Denote the projection (u,v) —> u by P. By the result of Corollary 13.21 

< C T . (3.14) 



E\u f 



L 2 (Q,T;H}{D t )) 



Write z e (t) as 

z e (t) = z e (0) - 

Then by (13.31) and (13.81) . when (h, 0) G X\ , we have the following estimate, for 
some positive constant C > independent of e 



A e z t (s)ds + / F e (s, x)rfs + / G t (s,x)dW(s). 



< 



< C 



a(Pz e (s), h)ds 
t 

\u e (s)\ H i (a 



ds 



(- P J A £ z e (s)ds + P J F e (s,x)ds,h 

e)dS + l \f( S ^ L2 W dS )\ h \ H o( D *> 
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Thus we have 

E -P [ A e z e (s)ds + P [ FJs,x)ds 2 <C T . (3.15) 

Jo JO W^(0,T;H-HD e )) 

Let M e (s,t) = J t s G e {s,x)dW{s). By Lemma 7.2 of PU and the Holder in- 
equality, we have 

E\PM e (s,t)\^ {De ) < E\PM e (s,t)\i>(D)<cE( f \ gi (r)\ 2 cQl dr) 2 

J s 2 

< K{t-s) f E\ gi (r)\ 4 cQl dr 

J s 2 

< K'(t-s) 2 

for t G [s, T], and for positive constants K and K' independent of e, s and t. 
Therefore 

E f \PM e (0,t)\t 2(Dt) dt<C T (3.16) 
Jo 



and for a G (|, |) 



; r /" T |PM e (0, t) — FM e (0, s) If 2/ n \ 
E/ / i )lL ^ dsdt<C T . (3.17) 



'o Jo r 

Combining the estimates fl3.14l) - fl3.17p with the Chebyshev inequality p^6| [19] . 
it is clear that for any 5 > there is a bounded set 

K s cX 

with X = L 2 (0,T;H^(D))f] (W /1,2 (0, T; H~ 1 (D)) + W a ' 4 (0, T; L 2 (D))) , such 
that 

PR G iQ} > 1 — 8. 
Moreover by the compact embedding 

L 2 (0, T; i? 1 ^)) D W h2 (0, T- H~ 1 {D)) C L 2 (0, T; L 2 (D)) n (7(0, T; H^^D)) 

and 

L 2 (0, T; i^p)) n H/ a ' 4 (0, T; L 2 (D)) C L 2 (0, T; L 2 (D)) n (7(0, T; H- l {D)), 

we conclude that ifj is compact in L 2 (0,T; L 2 (D)) n C(0,T; H~ 1 (D)). Thus 
{£(« e )} e is tight in L 2 (0, T; L 2 {D)) n (7(0, T; H-\D)). 

The proof is complete. □ 

Remark 3.4. When f = f(t,x,u e ) is nonlinear (i.e., it depends on u e ) but is 
also globally Lipschitz in u e , the results in Theorem \3.1\ and Corollary \3.2\ still 
hold. For example, see [11] for such SPDEs with stochastic dynamical boundary 
conditions. Moreover, by the Lipschitz property, we have \f(t,x,u^)\i?(p) < 
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Ct- Hence a similar analysis as in the proof of Theorem \3. 3\ yields the tightness 
of the distribution for u e in this globally Lipschitz nonlinear case. This fact 
will be used in the beginning of JSI to get the homogenized effective model when 
f = f(t,x,u t ) is globally Lipschitz nonlinear. 

In fact, in §02, we will also derive homogenized effective models for three 
types of nonlinearities f = f(t,x,u e ) that are not globally Lipschitz in u t . 

4. Two-scale convergence and some preliminary results 
In this section we present some basic results about the two-scale convergence 

mm- 

In the following we denote by C™(Y) the space of infinitely differentiable 
functions in R n that are periodic in Y. We also denote by L 2 er (Y) or Hp er (Y) 
the completion of C™(Y), in the usual norm of L 2 (Y) or respectively. 
We also introduce the space i7i, r (Y) /M, which is the space of the equivalent 
classes of u G ifi, r (F) under the following equivalent relation 

u^v^u — v = constant. 

Definition 4.1. A sequence of functions u e (t,x) in L 2 (Dj-) is said to be two- 
scale convergent to a limit u(t,x,y) G L 2 (DtxY) , if for any function ip(x,y) G 
C?(D T ,C™(Y)), 

lim / u e (t, x)ip(tj x, —)dxdt = — — j j u (t, x, y)tp{t, x, y)dydxdt. 



— I — c \ - 7 — / r \ ~ i 7 / t j* 

^°Jdt 6 \Y\JD T JY 

2—s 

This two-scale convergence is written as u e — > u. 

The following result ensures the existence of two-scale limit and for the proof 
see[2j[l2]. 

Lemma 4.2. Let u e be a bounded sequence in L 2 (Dt). Then there exist a 
function u G L 2 (Dt x Y) and a subsequence u ek with — ► as k — » oo such 
that u ek two-scale converges to u. 

Remark 4.3. Taking ip independent of y in the definition of two-scale conver- 
gence, then u t » u implies that u t weakly converges to its spatial average: 

u e (t,x) ->> u(t,x) = J u(t,x,y)dy. 

So, we see that, for a given bounded sequence L 2 (Dt), the two-scale limit 
u(t,x,y) contains more information than the weak limit u(t,x): u gives some 
knowledge on the periodic oscillations of u e , while u is just the average with 
respect to y. Another advantage of the usage of two-scale convergence is that 
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we do not need an extending operator such as in [331 03] in the homogenization 
procedure. For more properties of two-scale convergence we refer to [2J. 

The following result is useful when considering two-scale convergence of the 
product of two convergent sequences, see [2J, [12] . 

Lemma 4.4. Let v e be a sequence in L 2 (Dt) that two-scale converges to a 
limit v(x,y) G L 2 (Dt x Y). Further assume that 

lim / \v e (t, x)\ 2 dxdt = - — , / / \v(t,x,y)\ 2 dydxdt. (4.1) 
e ^°JD T \y\ Jd t Jy 

Then, for any sequence u e G L 2 (Dt), which two-scale converges to a limit 

u G L 2 (Dt xY), we have the weak convergence of the product u e v e : 

u e Ve j-Tj- J u{-,-,y)v(-,-,y)dy, as e -> in L 2 {D T ). 

Remark 4.5. Condition |^~Zp always holds for a sequence of functions (p(t, x, x/e), 
with (p(t,x,y) G L 2 (Dt', C per (Y)). Such functions v e are called admissible test 
functions. With the additional condition ( KT7p ; the two-scale convergence of v e 
is also called strong two-scale convergence [2j. 

Let u e be a sequence of functions defined on [0, T] x D e which is bounded in 
L 2 (0,T; H^(D e )). Then we have the following result concerning the two-scale 

limit of the bounded sequences u e and V x u e ; for the proof see [2j. 

Lemma 4.6. There exist u(t,x) G Hq(D t ), ui(t,x,y) G L 2 (D T ]H^ er {Y)) and 
a subsequence u ek with tk — > as k —>■ oo, such that 

u £k (t, x) x(y)u(t, x), k -> oo 

and 

V x M £fc x{y)[V x u{t,x) + V y ui{t,x,y)\, k -> oo 

where x{v) ^ s ^ e indicator function of Y* (which takes value 1 on Y* and 
value onY \ Y*). 



Since we consider the dynamical boundary condition, the technique of trans- 
forming the surface integrals into the volume integrals is useful in our approach. 
For this we follow the method of [55] (see also |14j ) for the nonhomogeneous 
Neumann boundary problem for an elliptic equation. 

For h G H^^i^dS) and F-periodic, define 

Ah = / K x ) dx - 

\y i Jas 



HOMOGENIZED DYNAMICS OF SPDES 



15 



Also define 

Thus, in particular Ai = jp-j and 



A = A, = ^, (4.2) 



where | • | denotes Lebesgue measure. 

For h G L 2 (dS) and F-periodic, define X e h G H~ 1 (D) as 



[ h{-)(p{x)dx, for <p G Hlip). 



Then we have the following result about the convergence of the integral on the 
boundary. 

Lemma 4.7. Let tp e be a sequence in Hq(D) such that ip e — ^ if weakly in 
Hl{D) ase^O. Then 



(X £ h ,(fe\D c ) -> A A / (fdx, as e—>0. 
Jd 

For the proof we refer to [14] , 

5. Homogenized macroscopic model 

In this section we derive the effective macroscopic model for the original 
model (12.11) . by the two-scale convergence approach. We first obtain a two- 
scale limiting model. Then the homogenized macroscopic model is obtained 
by exploiting the relation between weak limit and the two-scale limit. 

By the proof of Theorem 13.31 for any 5 > there is a bounded closed set 
K$ C X which is compact in L 2 (0, T; L 2 (D)) such that 

PR G K s } > 1 - 5. 

Then the Prohorov theorem and the Skorohod embedding theorem ([16]) as- 
sure that for any sequence {tj}j with e 3 - — > as j — > oo, there exist subsequence 
{ej(fc)}, random variables {w* } C L 2 (0, T; H £j{k) ) and u* G L 2 (0,T;H) de- 
fined on a new probability space (O* , JF* , P*), such that 



and 



£(m* ) = £({L. " 



m* — > u* m L (0, T; if) as k — ► oo, 



for almost all u; G fi. Moreover «,*. solves system (I2.ip - (I2.4I) with W replaced 



by Wiener process H 7 ^ defined on probability space (f2*,jF*,P*) with same 
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distribution as W. In the following, we will determine the limiting equation 
(homogenized effective equation) that u* satisfies and the limiting equation is 
independent of e. After this is done we see that C(u e ) weakly converges to 
C(u*) as e | 0. 

For u t in set Kg, by Lemma H~6l there is u(t, x) E Hq(Dt) and ui(t,x,y) E 
I?{Dt'i Hp er (Y)) such that 



and 



u e At,x) ^-4 x(y)u(t,x) 



- — — 2— s 

V x m . — > x(z/)[V x w(t,x) + Vj,tti(t,a;,j/)]. 



Then by Remark [O] 

M e . (t , x) — ^ ■7—j ■ / x(y) u (t, x)dy = $u(t, x), weakly in L 2 (Dt). 
\ Y \ Jy 

In fact by the compactness of Kg, the above convergence is strong in L 2 (Dt)- 
In the following, we will determine the limiting equation, which is a two-scale 
system that u and Ui satisfy. Then the limiting equation (homogenized effec- 
tive equation) that u satisfies can be easily obtained by the relation between 
weak limit and the two-scale limit. 

Define a new probability space (Qg, J-'g^g) as 

Qg = {cu E : u e (uj) E Kg}, 
Fs = {Fnn 5 :FEJ r } 



and 



Denote by Eig the expectation operator with respect to Pj. Now we restrict 
the system on the probability space (Qg,J-'g,¥'g). 

Replace the test function if in (13.71) by ip e (t, x) = <f>(t, x) + e<E>(i, x, x/e) with 
(j)(t,x) E C™(D T ) and &{t,x,y) E C™{D T ;C™ r {Y)). We will consider the 
terms in (13.71) respectively . 

By the choice of <p e and noticing that xd c weakly* in L°°(D), we have 

f(t,x)<f e (t,x)dxdt— / / XD e f{t,x)(p e (t,x)dxdt 



'0 JD e JO JD 

$ I I f{t,x)<f>(t,x)dxdt, e^O. (5.1) 

JO JD 
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And by the condition ( 12.51) 

ngi(t,x)tp e (t,x)dxdWi(t) = / / XD e 9i(t,x)(f e {t,x)dxdW 1 {t) 
_) e Jo Jd 



-»• / / g 1 (t,x)(j)(t,x)dxdW 1 (t), e^O m (5.2) 
Jo Jd 

Integrating by parts and noticing that u e converges strongly to i}u(t,x) in 

/. 2 UV). 




u e (t, x)tp € (t, x)dxdt = — I l u e (t, x)(p £ (t, x)dxdt 
'o J D t Jo J D e 

— — I / u e (t,x)cf)(t,x)dxdt — e / / u e (t, x)<b(t, x, — )dxdt 
Jo Jd Jo Jd e 

— > — / / l du(t,x)(f)(t,x)dxdt — I I 0u{t,x)4>{t,x)dxdt. (5.3) 
Jo Jd Jo Jd 

By the choice of ip e , 

V x (f>{t, x) + V„$(t, s, ~) ^ V x 0(t, x) + V„$(*, x, y), e -> 

and 

x 1 f 2 

]im\\V x (j)(t,x)+Vy$(t,x, -)||[i,2( DT )]n = —7 / V a 0(t, a;)+V y $(t, x, y) dydxdt. 

e \ Y I Jd t xY 

Hence by Theorem 14.41 we have 

nVu e (t,x)Vip e {t,x)dxdt = I f Vu e (t,x)(V x (f)(t,x) +V y ${t,x,-))dxdt 

= / [vu t (t,x)(y x <j) + Vy®(t,x,-))dxdt 
Jo Jd e 

-> J D f Y X(y)[Vxu(x) + V„ui(a;,3/)] [V^^x) + V v $(t,ar,y)]d2/da;tft 

1 ' T 



W\ I I I \yM x )+V v u 1 (x,y)][V x <j>(t,x)+Vy^t,x,y)]dydxdt. (5.4) 
I Jo JdJv 
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Now we consider the integrals on the boundary. First for a fixed T > 0, it is 
easy to see that 

/ u e (t,x)(p t (t,x)dxdt 

JdSe 



—e 



as, Jo 



u e (t, x)ip e (t, x)dtdx 
e(^X{, I u e (t,x)(p e (t,x)dt 



0, e^O. 



(5.5) 



And then 



eb 



u t (t, x)(p t (t, x)dxdt 



JdSe 



X\, / u e (t, x)(f e (t, x)dt 



D, 



MX 



(5.6) 



u(t, x)(fi(t, x)dxdt, e — > 0. 

'o JD 

By the same method as above and the condition (12.51) we have the limit of the 
stochastic integral on the boundary 

/ g 2 (t,x)ip t (t,x)dxdW 2 (t) 
o Jas e 

X I [ g 2 (t,x)<p(t,x)dxdW 2 {t), e -> 0, in L 2 (tt). (5.7) 
Jo Jd 

Combining the above analysis in (15. ip - (15. 71) and by the density argument we 
have 

d I I u(t, x)(f)(t, x)dxdt 



o Jd 



Kn I II [^xu{x)+V y u 1 {x,y)][V x (f){t,x)+V y ^{t,x,y)]dxdt 
\ Y \ Jo Jd Jy* 



-MX 



u(t, x)(fi(t, x)dxdt + $ 



'o Jd 

+W g 1 (t,x)4>(t,x)dxdWi(t) + X 
Jo Jd 



o Jd 

t 




f(t, x)(p{t, x)dxdt 
g 2 {t,x)4>(t,x)dxdW 2 (t) (5.1 



o JD 



for any G Hq(D t ) and $ G L 2 (D T ; Hp er (Y)/W). Integrating by parts, we see 
that (15. 8p is the variational problem of the following two-scale homogenized 
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system 

#du=[- div x A(V x u) - b<&\ x u + &f]dt + ■dgxdW^t) + \g 2 dW 2 {t), (5.9) 

[V x u + V«i] • v = 0, on - <9F (5.10) 
where i/ is the unit exterior norm vector on BY* — BY and 

M^xu) = ^J Nxu{t, x) + VyUtit, x, y)}dy, (5.11) 

with u\ satisfying the following integral equation 

J [V x u + V y Ui]Vy<bdy = 0, Mi is Y - periodic, (5.12) 

for any $ G Hq(Dt', H^ er (Y)). The problem f 1 5 . 1 2 j) has a unique solution for 
any fixed u, and so A(V x u) is well-defined. Furthermore A(V x u) satisfies 

and 

|^(0,OL W2 (D)|</3||^||i 2(D) (5.14) 

with some a, /3 > and any £, £x, £ 2 £ Hq(D). For more detailed properties of 
A(Vu) and (I5.12p we refer to [21]. Then by the classical theory of the SPDEs, 
HB|, fl5T9l- fl5TT0j) is well-posed. 
In fact A(Vu) can be transformed to the classical homogenized matrix by 

Ui(t, x,y)=^T t My) - e iV)> ( 5 - 15 ) 



where {ej}™ =1 is the canonical basis of ~R n and tOj is the solution of the following 
cell problem (problem defined on the spatial elementary cell) 

A y Wi(y) = in Y* (5.16) 

Wi — e^y is Y — -periodic (5-17) 
f)w ■ 

— - = on dS. (5.18) 
av 

Then a simple calculation yields 

A(Vu) = A*Vu 

with A* = (A*j) being the classical homogenized matrix defined as 

A$ = ± J WiMwjMdy. (5.19) 
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Then the above two-scale system (15.91) is equivalent to the following homoge- 
nized system, 

#du=[- div x (A*V x u) - bdXu + $f]dt + &gidWi(t) + Xg 2 dW 2 (t). (5.20) 

Let U(t,x) = l du(t,x). We thus have the limiting homogenized equation 

dU = [- d~ l div x {A*V x U)-bXU + df]dt + dg l dW l {t) + Xg 2 dW 2 (t). (5.21) 

And then it*, we have mentioned at the beginning of this section, satisfies 
(I5.2ip with W = (Wi,W 2 ) replaced by a Wiener process W* with the same 
distribution as W. By the classical existence result [16], the homogenized 
model (15.211) is well-posed. We formulate the main result of this section as 
follows. 

Theorem 5.1. (Homogenized macroscopic model) 

Assume that 112.5]) holds. Let u e be the solution of ll2.1}) - (2~4\ ). Then for any 
fixed T > 0, the distribution C(u e ) converges weakly to fi in L 2 (0,T; H) as 
e I 0, with \x being the distribution of U, which is the solution of the following 
homogenized effective equation 

dU = [-$- 1 div x (A*V x U) -bXU + df]dt + dg 1 dW 1 (t) + Xg 2 dW 2 (t), (5.22) 

with the boundary condition U = on dD, the initial condition U(0) = u^/d 
and the effective matrix A* = (A*j) being determined by A5.19\) . Moreover, the 

constant coefficient $ = is defined in the beginning of §2 and X = ^pp is 
defined in fiJJty . 

Proof. Noticing the arbitrariness of the choice of S, this is direct result of the 
analysis of the first part in this section by the Skorohod theorem and the 
L 2 (Qs)— convergence of u e on (Qg, J~s, P5). □ 

Remark 5.2. It is interesting to note the following fact. Even when the orig- 
inal microscopic model equation \2. 1\) is a deterministic PDE (i.e., g\ = 0), 
the homogenized macroscopic model 115.22}) is still a stochastic PDE, due to 
the impact of random dynamical interactions on the boundary of small scale 
heterogeneities. 

Remark 5.3. For the macroscopic system l}5. 22}) . we see that the fast scale 
random fluctuations on the boundary is recognized or quantified in the homog- 
enized equation, through the Hig 2 dW 2 (t) term. The effect of random boundary 
evolution is thus felt by the homogenized system on the whole domain. 
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6. Homogenized macroscopic dynamics for nonlinear 
microscopic systems 

In this section, we derive homogenized macroscopic model for the micro- 
scopic system (I2.ip - fl2.4l) . when the nonlinearity / is either globally Lipschitz, 
or non-globally Lipschitz. 

As Remark 13.41 has pointed out that if / is a globally Lipschitz nonlinear 
function of u e all the estimates in § [3] hold. In fact, similar results in §[5] on 
homogenized model also hold. In fact for / satisfying f(t, x, 0) = and 

|/(t,x,ui) - f(t,x,u 2 )\ < L\ui - u 2 \ 

for any t G R, x G D and Ui, u 2 G R with some positive constant L. Since u e — > 
"&u strongly in L 2 (0, T; L 2 (D)) and by the Lipschitz property of f(t, x, u) with 
respect to u, f(t, x, u e (t, x)) — » fit, x, u(t, x)) strongly in L 2 (0, T; L 2 (D)). (15.1 j) 
still hold for f(t,x,u e ). Then we can obtain the same effective macroscopic 
system as (15.221) with nonlinearity / = f(t,x,U): 

dU = [- $- l div x (A*V X U) - bXU + #f(t, x, U)] dt 

+dg 1 dW 1 {t) + \g 2 dW 2 {t). (6.1) 

For the rest of this section, we consider three types of nonlinear systems 
with / being non-global-Lipschitz nonlinear function in u e . The difficulty is 
at passing the limit e — > in the nonlinear term. These three types of nonlin- 
earity include: Polynomial nonlinearity; nonlinear term that is sublinear; and 
nonlinearity that contains a gradient term Vii e . We look at these nonlinearities 
case by case, and only highlight the difference with the analysis in §5. 

Case 1: Polynomial nonlinearity 

First we suppose / is in the following form 

f(t, x, u) = — a(t, x)\u\ p u (6.2) 

with < ao < a(t, x) < a\ for t G [0, oo), x G D. And p satisfies the following 
condition 

2 

p < , if n > 3; p G R, if n = 2. (6.3) 

n — 2 

For this case we need the following Weak convergence lemma from Lions [33J . 

Let Q be a bounded region in R x R n . For any given functions g e and g in 
V{Q) (Kp<oc), if 

\9e\LP(Q) < C> 9e — ► 9 in Q almost everywhere 

for some positive constant C, then g e — ^ g weakly in L P (Q). 
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Noticing that F e (t,x,z e ) = (f(t,x,u e ),0) and (F e (t, x, z e ), z t ) x o < 0, the 
results in Theorem 13.11 can be obtained by the same method as in the proof 
of Theorem 13.11 Moreover by the assumption (I6.3p . \f(t, x, w e )U 2 (D T ) < CV, 
which by the analysis of Theorem 13.31 yields the tightness of the distribution 
of u e . 

Now we pass the limit e — > in f(t, x, u e ). In fact, noticing that u e converges 
strongly to du in L 2 (0, T; L 2 (D)), by the above weak convergence lemma with 
g e = f(t, x, u e ) and p = 2, f(t, x, u e ) converges weakly to f(t, x, flu) in L 2 (D T ). 
Threfore by the analysis for linear system in §5J we have the following result. 

Theorem 6.1. Assume that Ii2.5\) holds. Let u e be the solution of li2.1\) - (2~4\ ) 

with nonlinear term f being A6'.2\) . Then for any fixed T > 0, the distribu- 
tion C(u e ) converges weakly to ji in L 2 (0,T; H) as e [ 0, with ji being the 
distribution of U , which is the solution of the following homogenized effective 
equation 

dU = [-V^div^VvU) -bXU + tif(t,x,U)]dt 

+$g 1 dW 1 (t) + \g 2 dW 2 (t), (6.4) 

with the boundary condition U = on dD, the initial condition U(0) = u$/d 
and the effective matrix A* = {A*j) being determined by Ii5.19\) . Moreover, the 

constant coefficient $ = is defined in the beginning of $2 and X = '-j^r is 
defined in \4.2\j . 

Case 2: Nonlinear term that is sublinear 

More generally, we consider / : [0,r]xDxR- > M. & measurable function 
which is continuous in (x, £) G D x R for almost all t G [0, T] and satisfies 

[f{t, x, 6) -f{t, x, &)] [6 -6] >0 (6.5) 

for t > 0, x G D and £x> 62 G R. Moreover, we assume that / is sublinear, 

\f{t,x,0\<g{t)(l + \Z\), £eR, t>0, (6.6) 

where g G Lf^JO, 00). Notice that under the assumption (16.51) and (16.61) . / may 
not be a Lipschitz function. 

By the assumption (16. 6p we can also have the tightness of the distributions 
of u e and also conclude that XD e f(t,x,u e ) two-scale converges to a function 
denoted by fo(t,x,y) G L 2 (Z>r xY). In the following we need to identity 
fo(t,x,y). 

Let G C^{D T ) and ^ G C%°{D T ; C™ r {Y)). And for re > let 

Ut,x) = <P(t,x) + KiPit,x,-). (6.7) 
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Then by the assumption (16.51) we have 
< 



with 



and 



n[f(t, x, u t ) - f(t, x, &)] [u e - Q dxdt 
/ x(-) [f(t,x,u £ ) -f(t,x,£ e )] [u e -Qdxdt 



= h — h,e — h,e ~ Ii,e + h,e 



X 

h,e = I x(-)f(^ x ^ e )u e dxdt 

D T e 



^ TvT / / fo(t,x,y)'&u(t,x)dydxdt, (6.8 
\ Y I io T Jy 

f X 

h,e = / x(-)f(t,x,u £ )£ e dxdt 



fo(t, x, y) [<j)(t, x) + Kip(t, x, y)] dydxdt, (6.9) 




D T JY 



h,e = / x(-)f{t,x,&)u e dxdt 



X(y)f(t, x, <j)(t, x) + Kip(t, x, y)) 




D T JY 



$u(t,x) dydxdt, (6.10) 



4,e = / x(-)f(t,X,£ e )Cedxdt 
JDt e 



1 

\y\ 



7T7T I I X{y)f(t,X,<j>(t,x) + K1p(t,X,y)) X 




' D T JY 

(t,x) + Ki/j(t,x,y)]dydxdt. (6-11) 

In (16.8l) - (l6.1ip we have used the fact of strong two-scale convergence of x(~) 
and f(t, x, £ e ), and the strong convergence of u e to flu. 
Now we have 



lim/ e 




D T JY 



[fo(t, x, y) - x(y)f(t, x,<p + AV>)] x) - <p{t, x) - m/>] dydxdt > 0, 
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for any G C$°(D T ) and i[> G C%°(D T ;C per (Y)). Letting -> dividing 
the above formula by k on both sides of the above formula and letting k — > 
yields 

/ / [fo(t,x,y)-x(y)f(t,x,^u)]^dydxdt<0 
Jd t Jy 

for any ^ G C£°(Dt; C per (Y)), which means 

fo(t,x,y) = x(y)f(t,x^u). 

Then by the similar analysis for linear systems in $5], we have the following 
homogenized model. 

Theorem 6.2. Assume that A2.5\) holds. Let u e be the solution of li2.1\) - (2~4\ ) 

with nonlinear term f satisfying A6.5\) and A6.6\) . Then for any fixed T > 0, 
the distribution C(u e ) converges weakly to \i in L 2 (0,T; H) as e [ 0, with \i 
being the distribution of U , which is the solution of the following homogenized 
effective equation 

dU =[- ®~ x div x (A*V X U) - bXU + 0/(t, x, U)] dt 

+$g l dW 1 (t) + Xg 2 dW 2 (t), (6.12) 

with the boundary condition U = on dD, the initial condition U(0) = u^/d 
and the effective matrix A* = (A*j) being determined by A5.19\) . Moreover, the 

constant coefficient ■d = is defined in the beginning of §2 and A = ^pp is 
defined in 

Case 3: Nonlinearity that contains a gradient term 

We next consider / in the following form containing a gradient term, 

f(t,x,u,'Vu) = h(t,x,u)-'Vu (6.13) 

where h(t, x, u) = (hi(t, x, u), ■ • • , h n {t, x, u)) and each hi : [0, T] x D x R — > K, 
i = 1, - ■ ■ , n, is continuous with respect to u and h(-, •, u(-, •)) G L 2 (0, T; L 2 (D)) 
for u G L 2 (0,T;L 2 (D)). Moreover assume that h satisfies 

(1) \(h(t,x,u) ■ Vu,v)l2\ < Co|VM|L2|t>| L 2 with some positive constant Cq. 

(2) \hi(t,x,£i) - hi(t,x,&)\ < k\£i - £ 2 | for £i, £ 2 e R, i = 1, • • • , n and k 
is a positive constant. 

Now we have 

\(F e (t,X,Z e ),Z e ) X o\ = \(h(t,X,U e ) ■ Vu e , Ue)L*\ < C \z e \ X o \z € \ X i ■ (6.U) 

By applying It 6 formula to |^ e |x°' we obtain 

d\z e (t)\ 2 x o+2(A e z e ,z e ) x odt = 2(F e (t,x,z e ),z e ) x odt + 2(G e (t,x)dW(t),z e ) x o + 

\G t {t,x)\ 2 Q dt. (6.15) 

*-9 
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By (EH, coercivity (E3D of a e 
(I6.15P with respect to t yields 



and the Cauchy inequality, integrating 



\ z t(t)\xO + a / \z t (s)\ x ids 



< \z \ x o + (2/5 + Ax (a)) / \z t {s)\ 2 x0 ds + 



(G e (s)dW{s),z e (s)) x o 



\G e (s)\ 2 cQ ds 



where A x is a positive constant depending on a. Then by the Gronwall lemma 
we see that (I3.9P and (13.101) hold. Moreover, the fact 

\h(t,x,u) ■ Vu\ L 2 < C \z e \ X }, 

together with the Holder inequality yields 

2 



E 



W 1 ' 2 {0,T;H~ 1 {D t )) 



(6.16) 



P [ A e z e (s)ds + P [ F e (s,x,z e )ds 
Jo Jo 

where P is defined in Theorem 13.31 Then by the same discussion of Theorem 
3.31 we have the tightness of the distributions of u e . 

Now we pass the limit e — > in the nonlinear term f(t, x, u t , Vw e ). In fact 
we restrict the system on (f^, JF^, P^). By the assumption (2) on h and the 
fact that u e strong converges to flu in L 2 (Dt), we have 

2 



lim 



dxdt = 0. 



h(t, x, u e (t, x)) — h(t, x, $u(t, x)) 
For any if) 6 C™{D T ), 

h(t, x, u e ) ■ Vu t ipdxdt 

h(t, x, u t ) — h(t, x, ■ Vu t ipdxdt + 

h(t, x, ■ Vu e ipdxdt 

— — / / h(t, x, flu) ■ x{y) \y xU + VyUi\ipdydxdt. (6.17) 
I * I Jd t Jy 

Combining with the analysis for linear system in §3 we have the following 
result. 

Theorem 6.3. Assume that Ii2.5\) holds. Let u t be the solution of ll2.1\) - (2~4\ ) 

with nonlinear term Ii6.13\) . Then for any fixed T > 0, the distribution C{u e ) 



D T 



2(> 



W. WANG & J. DUAN 



converges weakly to /i in L 2 (0,T; H) as e 10, with fi being the distribution of 
U = du which satisfies the following homogenized effective equation 

dU =[- $- l div x (A*V X U) - bXU + f*(t, x, U, V X U)] dt 

+$g 1 dW l (t) + Xg 2 dW 2 (t), (6.18) 

where the boundary condition U = on 3D, the initial condition £7(0) = u^/d, 
the effective matrix A* = (A*A is determined by $5.191) and f* is the following 
spatial average 

f*(t, x, U, V X U) = ±- J h(t,x,U) ■ x (y) [V-^xU + V„«i] dy 

with U\ being given by Ii5.15\) and x(v) the indicator function ofY*. Moreover, 
the constant coefficient d = is defined in the beginning of §2 and X = 
is defined in $4.2\ ). 

Remark 6.4. All the results in this paper hold when A is replaced by a more 
general strong elliptic operator d.iv(A e V 'u) , where A e is Y —periodic and satis- 
fies the strong ellipticity condition. 
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